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The anomalous scaling phenomena of three-dimensional passive scalar turbulence are studied 
using high resolution direct numerical simulation. The inertial range scaling exponents of the passive 
scalar increment and the scalar dissipation are obtained. The connection between the intermittency 
structure and the scaling exponent is examined and the structure instability of the high amplitude 
scalar dissipation is used to clarify previous experimental results for the scaling exponents. 
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The advection of a passive scalar field by a turbulent 
(or stochastic) incompressible velocity field is a process 
that exhibits cascade to small scales. The dynamics gov- 
erning a passive scalar is described by the following equa- 
tion, 



dT/dt - 



VT = kV'^T + f . 



(1) 



Here T is the passive scalar, k is the kinematic diffusiv- 
ity, f is a random forcing, and the advective velocity, u, 
is governed by the Navier-Stokes (NS) equations in three 
dimensions. 

The study of the dynamics in the passive scalar system 
has been one of the most active research areas in the field 
of fluid turbulence for the last decade In particular, 
based on a linear ansatz for the dissipation term, Kraich- 
nan Q obtains an explicit prediction of the anomalous 
scaling exponents for a class of passive scalar advected 
by a white and Gaussian incompressible velocity. The 
research has inspired a large number of papers |^ over 
the past two years, including numerical simulations ^] 
and analysis of experiments pj . 

The passive scalar convected by Navier-Stokes turbu- 
lence, however, is a more difficult problem where the ve- 
locity is far from white and Gaussian. No inertial range 
scaling based on the passive scalar dynamical equation 
(1) has been developed. The local similarity theory for 
the passive scalar, an extension to Kolmogorov 1941 sim- 
ilarity theory for the velocity field was studied by 
Obukhov 1^] and Corrsin ||] . For fluid flows at very high 
Reynolds number and with Prandtl number, Pr — i^/k, 
close to unity, Obukhov's theory assumes the existence 
of universal statistics of fluctuations at so-called inertial- 
range scales L ^ r ^ 77, where is the kinematic vis- 
cosity, L and rj are the characteristic length scales for 
the large scale and the dissipation scale, respectively. 
This hypothesis has a series of implications, including 
the 2/3 law for the passive scalar fluctuations in the iner- 
tial range: (AT^) ~ (e)'^^^ {N)r^/^ , where TV = k(VT)2 
is the scalar dissipation function, e is the velocity dis- 
sipation function, AT^ — T{x -I- r) — T{x) is the scalar 
increment and (•) denotes an ensemble average. In gen- 



eral, if there is a scaling range, the pth order structure 
function, Sp{r) = (AT^), should have a scaling relation: 
Sp{r) ~ r^p, with Zp defined as the pth order scahng 
exponent. 

Taking account of the intermittency correction J9] and 
using the refined similarity hypothesis (RSH) ||l^,|ll | , the 
phenomenological scaling models for the passive scalar, 
including the bi-variate log-normal model by Van Atta 
Jl2t , the (3 model jl^ by Frisch et. al. and the bi-variate 
multifractal model by Meneveau et. al. |lj] , lead to ana- 
lytical predictions of anomalous scaling exponents. In a 
recent note, we have developed a bi-variate log-Poisson 
model p^ ; the resulting scaling exponents agree well 
with experimental measurements. 

The challenge in the real-life experimental measure- 
ment is to minimize the interference of two probes which 
measure velocity and scalar quantities simultaneously at 
the same position |lj] . Our research is motivated by the 
paper by Sreenivasan and Antonia p6[ | who pointed out 
that two existing experimental measurements |17 ij] do 
not provide a convergent scaling relation for the passive 
scalar system. In particular, the scaling exponents in |lj] 
were obtained using the measured dissipation exponents 
by invoking the RSH, showing a saturation for p > 6. 
The scaling exponents from another experimental mea- 
surement p8[ fall between above two results. 

Direct numerical simulation of the passive scalar equa- 
tion (1) and the Navier-Stokes equations were carried out 
simultaneously. A mesh size of 512'^ was used in a cyclic 
cubic box for homogeneous isotropic turbulence. In or- 
der to maintain statistical steady states, both the veloc- 
ity field and the passive scalar field were forced for /c < 3 
||l9|,0. The forcing scheme keeps the total energy in 
the first two shells (1 < fc < 2 and 2 < k < 3) consistent 
with k~^/^ . In addition, Fourier modes in each shell have 
equal energy and the phase of each mode is randomized. 
For simplicity the Prandtl number is fixed to be unity for 
this study. The analysis was carried out for the statisti- 
cal steady states at Taylor microscale Reynolds number 
Tlx = 220. A spatial averaging over the whole physical 
domain and a time average over approximately 10 large- 
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eddy turnover time were used to replace the ensemble 
average. 
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FIG. 1. Numerical verification of Kolmogorov's 4/5 law 
and Yaglom's 4/3 law. The solid line is for Kolmogorov's 
4/5 law and the dashed line is for Yaglom's 4/3 law. The 
signs are from DNS data (A for — ((Aitr)^)/((e)r) and • for 
-{{Aur{ATrf)/{{N)r). Here = {u^/{e)y^'^ is the Kol- 
mogorov dissipation length, In the inset shows the scalar 
spectrum compared with Oboukov's —5/3 law for the iner- 
tial range scaling. 

At very high Reynolds number, Kolmogorov's 4/5 law 
pl[ is exact in the inertial range: {{Aur)^) = {{u{x + 
r) — u{x))^) = — 4/5(e)r. Similarly there is also an ex- 
act scaling relation for the third order cross correlation 
given by Yaglom ||: {Aur{ATr)'^) = -4/3(7V)r. To 
demonstrate the quality of the inertial range scaling of 
our DNS data, in Fig. 1 we present —{{Aur)^)/r{e) and 
— {Aur{ATr)'^)/r{N) as functions of r. It is seen that 
both quantities show a narrow inertial scaling range (less 
than one decade and they are also displaced with respect 
to each other). On the other hand, although DNS was 
carried out at moderate Reynolds numbers, the values of 
the dimensionless constants agree well with theoretical 
predictions (4/5 and 4/3 for velocity and passive scalar, 
respectively). It was reported in reference Q that these 
dimensionless constants were smaller in the experimental 
measurements (for a circular jet with TZx — 220 and the 
atmospheric surface layer with TZx = 7200) than the the- 
oretical values. In the inset, we show the passive scalar 
spectrum, Ep(k), as a function of the wave number, k. 
The inertial range with a slope slightly less than —5/3 
could be identified. This smaller than —5/3 exponent 
was previously observed experimentally [Q. A bump 
between the inertial and dissipation ranges perhaps is 
associated with the "bottleneck effect" psj . 

In Fig. 2, we plot structure functions, (lAT^I^), as 
functions of r for p = 2, 4, 6 and 8. The power law iner- 
tial range can be identified for 0.2 < r < 1. From data 
analysis we notice that in fact the passive scalar incre- 
ment displays a better scaling relation than the velocity 
increment , and therefore it is not necessary to invoke 



the Extended-Self-Similarity (ESS) technique (2^. As a 
matter of fact we have tried ESS by plotting {\ATrY') as a 
fimction of {Aur(ATrY) . The error bar for the extracted 
scaling exponents in the inertial range seems bigger than 
the result using the original structure functions. In inset 
(a), we show the local slope, Zp{r) — d\og Sp{r) / dlogr , 
as a function of r for the same p in the structure function. 
The local slope was calculated using a least square fit to 
a power law for every 3 neighboring points. A flat region 
for each p can be seen for Zp{r), supporting the existence 
of inertial scaling range. In inset (b), we present the two 
point correlation function of the passive scalar dissipa- 
tion function, {N{x)N{x + r)), as a function of r. Again 
a scaling relation for the dissipation correlation is found: 
{N{x)N{x + r)) - r-f" with ^ « 0.25. The intermittency 
parameter fi agrees well with previous experimental re- 
sults 
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FIG. 2. The structure function Sp{r) as a function of r. In 
the inset (a) shows the local slope Zp{r) as a function of r and 
in the inset (b), we present the two-point correlation function 
of the passive scalar dissipation. 

In Fig. 3, we show the normalized probability den- 
sity function (PDF), {D)P{D), as a function of the nor- 
malized fluctuation dissipation: D/{D), where D is the 
dissipation function. It is clear that both the velocity 
dissipation and the scalar dissipation fields are strongly 
intermittent |2^ . We have also compared the PDFs with 
the log-normal distribution by looking at the statistics of 
WD = {InD - {InD)) / {{InD - {lnD})^y/'\ It is found 
that both fields are quite close to the log-normal val- 
ues for low order moments. From the plot, we also note 
that the PDF of the scalar dissipation has a wider tail 
than that of the velocity dissipation, indicating that spa- 
tially there are more large amplitude events in the passive 
scalar dissipation than the velocity dissipation. In the in- 
set we show the flatness of the passive scalar increment 
and the longitudinal velocity increment as functions of 
the separation. It is seen that the flatness of the scalar 
increment is larger than that of the longitudinal veloc- 
ity increment, except when r is in the large scale region 
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where both fields are essentially Gaussian and their flat- 
nesses approach the Gaussian value of 3. In particular, 
when r is in the dissipation range, the flatness of the pas- 
sive scalar increment is significantly larger than that of 
the velocity increment, in agreement with the PDFs. 
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FIG. 3. Normalized PDF, {D)P{D) as a function of tlie nor- 
malized dissipation. Here D is the dissipation function (the 
solid line is for the scalar dissipation and the dotted-dash line 
for the velocity dissipation). The dotted line is for the scalar 
dissipation function with Gaussian statistics. In the inset we 
present the flatness of the scalar increment (□) and the ve- 
locity increment (x) as a function of r (lattice unit). 




FIG. 4. The scaling exponent, Zp, as a function of p. The 
• sign is for DNS and the □ for the experiment data in [ p^ . 
The solid line is for a bi-variate log-Poisson model in In 
the inset shows the scaling exponents at two different time 
steps without using time averaging. The data from are 
represented by the o sign. 



The scaling exponent, Zp, as a function of the order 
index, p, is shown in Fig. 4. The • sign is for the result 
from DNS measurement and the □ sign for the result 
from the experiment by Antonia et al. The solid 

line is from a bi-variate log-Poisson model in jl^ assum- 
ing that the correlation coefficient between velocity field 



and the passive scalar field is zero. The DNS scaling ex- 
ponents are obtained using the flat region in the the local 
slope plot as shown in the inset (a) of Fig. 2. Some typ- 
ical scaling exponents and the corresponding errors are 
listed here: zo.5 = 0.165 ± 0.007, zi = 0.354 ± 0.012, Z2 = 
0.606 ± 0.019, Z3 = 0.794 ± 0.025, = 0.943 ± 0.030, = 
1.180±0.041,Z8 = 1.376±0.050 and zw = 1.546±0.063. 
It is noticed that the second order scaling exponent devi- 
ates from Obukhov's 2/3 law by 15%. This result agrees 
well with previous experimental results ||l^,|lJ,Q. The 
fourth order scaling exponent in DNS also agrees quite 
well with the theoretical prediction: 24 — 2z2 — /i in [Q, 
a result primarily for the white and Gaussian velocity. 
In general, all scaling exponents from DNS coincide well 
with the phenomenological theory and the real-life 
experiment [ pT| . 

We point out that the DNS results shown in Fig. 4 were 
obtained using a spatial averaging over the whole simu- 
lation domain (512'^) and a long-time averaging (~ 10 
large-eddy turnover times). In fact, the scaling exponent 
for each single time frame displays intense fluctuation 
with time, much larger than the fluctuation of the expo- 
nents for the velocity increment [ p^ using a spatial aver- 
aging of the same size. To demonstrate this fluctuation, 
in the inset of Fig. 4, two typical scaling exponents with- 
out time averaging are shown. The time difference be- 
tween two frames is about 0.2 large-eddy turnover time. 
Although the curves for p < 4 agree quite well, they 
arc qualitatively different for larger p. A large variation 
of the corresponding PDFs with time has also been ob- 
served. One set of the scaling exponents (□) clearly es- 
tablishes a saturation when p > 7, indicating a possible 
upper bound for the most intensive events. The satura- 
tion of scaling exponents for the passive scalar has been 
observed in (shown also in the inset by the o sign) 
and . The scaling exponents in another group ( x ) are 
larger than the time- averaged values in Fig. 4 for p > 5. 
Such strong fluctuations could be a reason for the satu- 
ration displayed in [Q where the results were obtained 
based on 1-D cuts through the held and may not capture 
the strongest events with a limited sample. On the other 
hand, the experimental scaling exponents were obtained 
through the RSH and it can not be ruled out that the 
saturation could be due to shortcomings of the RSH at 
high values of p. 

To understand the physics behind the large time fluctu- 
ation of the scaling exponent, we have studied the inter- 
mittency structures of the passive scalar fleld, in particu- 
lar the large amplitude events of the dissipation function 
and the scalar derivative. The latter is directly associ- 
ated with the scalar incr ement: AT^ = d{T/dx)dx. 
Using visualization technique, we flnd that with increas- 
ing of amplitude, the iso-surface of the scalar dissipation 
changes from fragment-like to sheet-like (shown in Fig. 
5). The dynamical evolution of the sheet-like structures, 
such as formation and annihilation, clearly plays a key 
role in the time dependence of the scaling exponents, sim- 
ply because the large amplitude events are the predomi- 
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nate effects on the high order scalar increments and their 
exponents. From the structure stabihty point of view, we 
suspect that the sheet-hke structure are more unstable 
to perturbations, such as forcing, than the filament-like 
structure [^0|, which is a characteristic vortex structure 
for large amplitude events in three-D turbulence. This 
agrees with the observation that the velocity scaling ex- 
ponent seems to be less time dependent. The instability 
of the sheet-like structures possibly causes the annihila- 
tion of the very high intermittency structures at a specific 
time, leading to a saturation of the scaling exponents. In 
addition, we know from Fig. 3 that spatially there are 
more high amplitude events in the scalar field than in the 
velocity field, indicating that the generation and the an- 
nihilation of the sheet-like structure is more frequent and 
causes intense fluctuations in the statistics of the scalar 
field. 




FIG. 5. Iso-surface of the scalar dissipation at the value 
A'' = 10{N). The resolution of the plot is 64^ and the data 
shown represent the neighborhood of the overall maximum of 
scalar dissipation. 

The current finding of the time dependence of the scal- 
ing exponents might connect with the idea of the non- 
universality of scaling exponents for the passive scalar 
system ||2^ in the sense that the detailed intermittency 
structure strongly affects the exponents. On the other 
hand, so far little is known about how the velocity 
stretching affects the dynamical evolution of the sheet- 
like structure in the scalar field. We feel that the fun- 
damental physics of the scaling dynamics for the passive 
scalar turbulence is still missing and a more detailed and 
careful study is much needed. 

We thank R. Kraichnan, C. Meneveau, M. Nclkin, Z.-S. 
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ical simulation was carried out at the Advanced Comput- 
ing Laboratory at Los Alamos National Laboratory. 
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